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Abstract 

In nuclear chiral perturbation theory (%PT), an operator is defined in a space with a cutoff which 
may be varied within a certain range. The operator runs as a result of the variation of the cutoff 
[renormalization group (RG) running]. In order for xPT to be useful, the operator should run in 
a way consistent with the counting rule; that is, the running of chiral counter terms have to be 
of natural size. We vary the cutoff using the Wilsonian renormalization group (WRG) equation, 
and examine this consistency. As an example, we study the s-wave pion production operator for 
NN — > dn, derived in xPT. We demonstrate that the WRG running does not generate any chiral- 
symmetry-violating (CSV) interaction, provided that we start with an operator which does not 
contain a CSV term. We analytically show how the counter terms are generated in the WRG 
running in case of the infinitesimal cutoff reduction. Based on the analytic result, we argue a range 
of the cutoff variation for which the running of the counter terms is of natural size. Then, we 
numerically confirm this. 
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I. INTRODUCTION 



In nuclear chiral perturbation theory (%PT), operators (e.g., nuclear force, electroweak 
current and pion production operator) are derived from a chiral Lagrangian following a 
counting rule. Many processes in few-nucleon system have been successfully described by 
these chiral operators, establishing the validity of this approach. The chiral nuclear 

operators are defined in a model space with a certain cutoff. The cutoff has a physical 
meaning and its choice is not arbitrary. Given a long-range (e.g., pion exchange) mechanism, 
the cutoff should be smaller than the scale where the details of shorter-range mechanisms are 
resolved. For example, in describing a nuclear operator with pion-exchange mechanisms plus 
contact terms, the cutoff should be smaller than the p meson mass. Also, the cutoff should 
not be so small that the long-range mechanism is not fully taken into account. Although not 
arbitrary, the cutoff still may be varied within a physically reasonable region. As a result 
of the variation of the cutoff, an operator (more specifically, couplings of counter terms) 
runs so that observables are cutoff-independent [renormalization group (RG) running]. A 
question here is whether the running of the counter terms is consistent with the counting 
rule. Stated differently, we wonder whether the counter terms run with keeping the size of 
the couplings natural [(9(1)]. This consistency is a necessary condition for the counting rule 
to be useful because, if not satisfied, the counting rule does not correctly reflect the ordering 
of the importance of the counter terms. 

In order to address the question concerning the consistency between the RG running 
and the counting rule, we propose to use a RG equation which controls the running of an 
operator so that Green's functions are unchanged. Wilsonian RG (WRG) equation is such 
an equation. For studying the RG running, we consider it the most appropriate to use the 
WRG equation because we have the following beneficial points: (i) The WRG equation is 
consistent with the derivation of the effective Lagrangian with a path integral, see Sec. [Til 
(ii) It is guaranteed that no chiral-symmetry-violating (CSV) interactions are generated 
by the WRG running, 1 provided that no CSV term exists (e.g., the operator vanishes 
at threshold in the chiral limit) in both a transition operator and a nuclear force before 
the RG running, see Sec. IIV( (iii) Since we have the RG running of the operator which 
correctly reflects the high- momentum states integrated out, we can study the convergence 
of the chiral expansion of the RG running. Our usage of the WRG equation is different 



from Refs. 



in which the scaling behavior of each interaction near the fixed point is 



identified with power counting. What we would like to do is also different from Refs 



. BB 



1 In Sec. [Ill we will explain in more detail what we mean by "no CSV terms are generated in the WRG 
running" . 
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where the renormalizability of the leading order (LO) chiral NN potential over a very wide 
range of the cutoff is examined. We apply the WRG equation to a set of operators derived 
with an established power counting, and then examine the internal consistency between the 
RG running and the counting rule over a physically reasonable range of the cutoff. In the 
above RG analyses^, Q, 0, 0], the authors are not primarily interested in the actual size of 
the running for a physically reasonable range of the cutoff. We are interested in the actual 
size of the RG running and its consistency with the counting rule. One may naively expect 
that this kind of calculation only generates the counter terms which scale as oc 1/A (A : 
cutoff). We will semi-analytically show that this expectation is not always the case, and find 
a size of the variation of A for which this expectation holds true. Also, we argue that the 
RG running is consistent with the counting rule if we vary the cutoff smaller than this size; 
otherwise, the RG running could be divergent. This argument is also consistent with the 
divergent RG running of some contact A^A" interactions which was found in Ref. j^J (e.g., 
see Fig. 2 of the paper). The physically reasonable range of the cutoff is consistent with the 
size of the variation of A which we semi-analytically find. 

We quantitatively confirm the above qualitative analysis by a numerical calculation. 2 In 
the numerical calculation, we use a phenomenological nuclear force rather than a xPT-based 
nuclear force. A reason is that chiral ATAT-potentials available^] have relatively small cutoffs 
(A ~ 500 MeV) and are not very appropriate to study the WRG running in the cutoff range 
considered in this work (500 MeV < A < 800 MeV). Phenomenological A^A^-potentials such 
as the CD-Bonn ATAf-potential [10] which we use are not based on the chiral Lagrangian. 
However, these potentials have been often used with chiral operators (electroweak currents, 
pion production operators) to calculate matrix elements (hybrid approach). This hybrid 
approach has been used extensively, through which its usefulness has been established 
Despite the phenomenological success of the hybrid approach, there is still a concern about 
the consistency between the nuclear force and the chiral transition operators. In this context, 
the present RG analysis may provide an interesting test for the hybrid approach. This is 
because if a CSV term is generated in the WRG running, then the result signals that the 
nuclear force contains a numerically significant CSV component and one has to doubt the 
validity of the hybrid approach. 

This paper is organized as follows. In Sec. [TTJ we first develop a formal apparatus to 
address the question of how to evolve the operator for changing A, using the WRG equation. 
In Sec. PH we present expressions for the chiral next-to-leading order (NLO) s-wave n- 



2 Similar analyses addressing the RG running and the naturalness of the counter terms in the NN interaction 
are found in Refs. 0, 0] ; the WRG equation is used in Ref. Q while a unitary transformation is used in 
Ref. 
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production operator 11] to which we apply the WRG analysis as a demonstration, and 
address the question of the consistency. The operator we use is based on the heavy-baryon 
%PT with a revised power counting scheme 12j which treats the large incoming nucleon 
momentum as a separate large energy /momentum scale (see, e.g., Ref. [lsj for a treatment 
using the original Weinberg counting [l4|). We use the revised power counting throughout 
this work unless otherwise specifically stated. In Sec. HVl we reduce the (sharp) cutoff of the 
chiral 7r-production operator by the infinitesimal amount using the WRG equation and the 
chiral LO nuclear force. We obtain analytic expressions which show that the WRG equation 
indeed induces the running of the chiral counter terms; no CSV terms are generated. In 
Sec. |Vl in order to quantitatively address the RG running, we reduce the cutoff of the ir- 
production operator from A = 800 MeV down to A = 700, 600 and 500 MeV using the 
WRG equation with a phenomenological nuclear force. As a result we obtain numerically 
the corresponding effective operators. Then, we try to reproduce the effective RG low- 
momentum operator 3 by calculating another operator consisting of the chiral NLO operator 
and the chiral counter terms; the chiral counter terms are introduced following the counting 
rule. The counter terms are accompanied by unknown coefficients, the so-called low-energy 
constants (LECs), which are determined by fitting to the low- momentum operator. We will 
show that the running of the operator with A is, to a very high precision, captured by the 
lowest order counter term. Finally, we summarize our findings in Sec. IVII Some calculations 
are relegated to an Appendix. 



II. WILSONIAN RENORMALIZATION GROUP EQUATION FOR TRANSI- 
TION OPERATOR 

The WRG equation is an equation which is designed to control the running of a set of 
operators, as a result of the cutoff reduction, so that Green's functions are unchanged. At 
first, we state why we use the WRG equation to reduce the cutoff. An effective Lagrangian 
can be obtained formally via a path integral formulation based on the Lagrangian of the 
underlying, more fundamental theory. One integrates out the high energy degrees of freedom 
using the path integral. When integrating out the high momentum states of the nucleon 
in the heavy-baryon xPT Lagrangian, we can also use the path integral. This procedure 
is equivalent to solving the WRG equation derived below; we will state more on this later. 
Thus one way of reducing the cutoff consistently with effective field theory is through using 
the WRG equation. It is not a priori guaranteed that the WRG running of the counter terms 

3 In this work, we will use the term "low-momentum interaction" to refer to an effective interaction obtained 
with the Wilsonian rcnormalization group equation. 
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is consistent with the counting rule, and thus we consider it worthwhile examining. We will 
see that the WRG running is consistent when the cutoff is changed within a physically 
reasonable range; otherwise, not necessarily consistent. 

Contrary to the approach in the above paragraphs, some previous works did not use a 
RG equation for studying RG running. Instead, for a given cutoff, the (leading) counter 
term was fixed so that some observables were reproduced. Repeating this procedure over a 
certain range of the cutoff leads to the RG running of the counter term. In this procedure, 
it is assumed that the running of the (leading) counter term simulates the high momentum 
states integrated out. However, one could fix the coupling of a higher-order counter term, 
instead of the leading term, so that some observables were reproduced. Without using a 
RG equation, it is impossible to know which counter term is running. Furthermore, our 
approach based on the WRG equation is more advantageous than the above RG analyses 
on the following two points: (i) it is guaranteed that no CSV operator is generated in the 
RG running, provided that no CSV term exists before the RG running; (ii) the RG running 
of the operator which exactly reflects the high momentum states integrated out is at hand, 
and thus we can examine whether the RG running is captured by a series of the chiral 
counter terms consistently with the counting rule. We detail the above statement (i) in the 
following. The chiral operators in the nuclear x?T are not chiral invariant in the sense that 
the original Lagrangian is invariant under chiral transformations. We take advantage of the 
chiral symmetry by deriving the operators from the chiral Lagrangian; the parameterization 
of the operators is given by the chiral Lagrangian. After the WRG running, if the operators 
are still accurately parameterized by the same parameterization, we state that no CSV 
interactions are generated. We do not address the chiral symmetry invariance of the original 
chiral Lagrangian in the RG running. Here we are interested in nuclear operators based on 
%PT, which is actually used in practical calculations, and their WRG running. As far as 
we examine the WRG running of the xPT-based operators with the cutoff, no CSV term is 
generated. 

As stated in the previous paragraph, the WRG equation can be derived using the path 
integral. The present author derived the WRG equation for the NN interaction in this way 
in Ref . 0] • The WRG equation for a transition operator (ir production operator in our case) 
can also be derived in essentially the same way, whereas we can also derive it in a simpler 
way as detailed in Ref. Il5] (Appendix A of that reference). 4 Here we derive the WRG 
equation following Ref. |l5| for simplicity without using the partial wave decomposition that 



In Ref. [3] , the WRG equation was derived as a sufficient condition of the cutoff independence of a matrix 
element. Recently, it was shown that the same WRG equation is also necessary-sufficient condition, if one 
imposes the cutoff independence of the five point Green function [l6| . 
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was used in Ref. [15(. We start with a matrix element in which the transition operator is 
defined in a model space spanned by plane wave states of the two-nucleon system. The 
maximum magnitude of the relative momentum in the model space is given by the cutoff, 
A. The matrix element of an operator O is given by 

(p'\o\ p ) = [ A dk f A dk' ^(fe') o{k\ k) i; p (k) , (i) 

Jo Jo 

where ip p {k) [ippi(k')} is the wave function for the initial (final) two-nucleon state. The wave 
functions are derived from a low-momentum AW interaction with the same cutoff A. The 
quantity p (p') is the on-shell relative momentum for the initial (final) two-nucleon state. 
The on-shell momentum is related to the energy (E) for the relative motion of the two 
nucleons through p = \p\ = V ME, where M is the nucleon mass. We denote the transition 
operator 0(k', k), where k (k') is the relative off-shell momentum of the two-nucleon system 
before (after) the interaction O. 

We differentiate both sides of Eq. ([1]) with respect to A and impose the renormalization 
condition that the matrix element is invariant under cutoff changes, i.e., d(0)/dA = 0. This 
gives the WRG equation for the low-momentum transition operator (with arguments now 
explicitly shown), 

dO A (k',k;p',p) 
dA 

M f dQ A ( OA(k',A;p',p)V A (A,k;p) V A (k', A; p')Q A (A, k;p',p) \ 

J (2tt) 3 { l-p 2 /A 2 l-p'1/A 2 )' U 

with A = A/A. The low-momentum AW-potential, V A (k',k;p), evolves according to the 
WRG equation for the NN potential . The low-momentum operator acquires a depen- 
dence on both the initial and final on-shell momenta, p and p', as indicated by Eq. (T5]). 

The WRG equation is solved in order to derive A from O = A with a cutoff A(> A). 
The solution (in integral form) of the WRG equation is (for later convenience now given 
with partial wave decomposition), 

of' a) = Jof a) + of a) — l —xv^ + v^ x l -^of a) 

+ vf>A KgrOf a) J-rrAvi^ n , (3) 

E> - Hf\ A E — XH^ j K J 

where and are the full Hamiltonian and the AW-interaction for partial wave a, 
defined in the model space with cutoff A. The operator O^f'^ generates a transition from a 
partial wave a to (3. The projection operators 77 and A are defined by 



V 



f k 2 dk 




k)(k 


1 2tt 2 



k < A , (4) 
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A 



f k 2 dk 




k)(k 


1 2vr 2 



A < k < A 



(5) 



where \k) represents the radial part of the free two-nucleon states with the relative mo- 
mentum k. Equation (J3J) is the same as for the effective operator in the Bio ch- Horowitz 
Similarities between the projection formalisms (e.g., Bloch-Horowitz and 
isms) and RG techniques have been explored previously in other contexts, 



formalism 
Lee-Suzuki forma 
see, e.g., Refs. [18 



19] . 



III. CHIRAL s-WAVE PION PRODUCTION OPERATOR 



Here we present the xPT-based s-wave pion production operator for the AW — > dn 
reaction near threshold. We will start our RG analysis with this operator, and examine 
the consistency between the RG running of this chiral operator and the counting rule. 
We use the operator from Ref. ll|, which was derived using the modified power counting 
rule proposed in Ref. [I2j. In this counting the large relative momentum (p = y/Mm n ) 
of the incoming nucleons is counted as an additional large energy-momentum scale of the 
problem, eading to an expansion in \ ~ y^ff rather than in the ^ of the original Weinberg 
counting Jjj]. The leading-order (LO) operators are given by the (nucleon recoil) one-body 
operator and rescattering via the Weinberg- Tomozawa (WT) interaction. At the next-to- 
leading order (NLO), several loop diagrams start to contribute, the sum of which does not 
vanish in the chiral limit when sandwiched between wave functions 20 L leading to a divergent 
matrix element. A solution to this dilemma was proposed in Ref. ll|, where it was shown 
that the WT term and its recoil correction, taken together with a pion exchange extracted 
from the initial or final state wave functions, contribute an irreducible NLO diagram that 
exactly cancels the offending divergence of the pion loop diagrams. 5 We are left with WT 
rescattering, with its energy dependence replaced by the on-shell value, and no CSV term. 
The result is that, up to NLO, we consider the (nucleon recoil) one-body plus a modified WT 
rescattering term (WT'), the latter a factor 4/3 stronger than the original WT term 11]. 
These operators are given (in momentum space) by 



O 



WT' 



o 



IB 



9_A^q abcbc 

-ig A w q C27i) 3 



ct x ■ (k' x - ki) 



<x 2 ■ (k' 2 - k 2 



ml 



-T, 



(k[ 



(k[ - fcO 2 m; 2 + (k' 2 
(k' 2 -k 2 )a l -(k 1 + k' 1 )) 
ki) <r 2 • (k 2 + k' 2 )) , 



k,)* 



T, 



(6) 



(7) 



In an irreducible diagram, the nucleons of any two-nucleon cut are off shell by ~ m ff , i.e., an irreducible 
diagram cannot be split into smaller diagrams with all external nucleons on-shell. 
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where m' n 2 = \m 2 K and u q is the energy of the emitted pion. The momentum for i-th nucleon 
in the initial (final) state is denoted by ki (feQ. The axial- vector coupling constant and the 
pion decay constant are denoted by qa and f w , respectively. We have employed the so-called 
fixed-kinematics approximation, in which the energy transfer is equally shared between the 
incoming nucleons and fixed to the threshold value m n /2, where m n is the pion mass. Other 
choices are possible 21[, but we will take this simple prescription here, and relegate an 
investigation regarding this issue to future work. 

When we let the one-body operator run according to the WRG equation, we obtain a 
low-momentum operator with a kink structure. (See Fig. [TJ) The origin of the kink is high 
momentum components of the bare one-body operator that are integrated out. We explain 
here this point more using the WRG equation [Eq. (j2j)]. When the cutoff is reduced by 5 A, 
the running of the operator due to the momentum shell of the one-body operator integrated 
out is given up to the order of 5 A as 

sn (V 1, > \ M f f OiBA(k',A)V A (A,k;p) V A (k',A;p')Q 1BA (A,k) \ 
50 A (k,k;p,p) = MJ — ^ ^ + j— ^ j 5 A ,(8) 

where the relative nucleon momenta before and after the insertion of O^b are k = {k\ — k^)/2 
and k' = (k[ — k' 2 )/2, respectively, and they are related to the pion momentum (q) through 
k' — k — |. The first term in r.h.s. of Eq. (JSl) is non-vanishing only when k' + 1 = A because 
of the 5-f unction in Eq. (J7|). This means that the shift of the operator [50 A (k', k;p',p)} is 
generated only for |fc'| > A — |||. Similarly, the second term in r.h.s. of Eq. (jBJ) induces 
the running of the operator only for |fe| > A — |||; the kink structure is created in this 
way. On the other hand, if Oib in Eq. (jSJ) is replaced by Owt' [Eq. ((SD], then the shift 
of the operator [50 A (k',k;p',p]] is induced for all values of k and k' because Owt' does 
not contain the 5-function. Thus the kink structure does not mean the strong dependence 
of the operator on the cutoff. Rather, it originates from the ^-function in the one-body 
operator, or in other words, from the way we have chosen to define the transition operators 
and wave functions. We could have chosen a set of transition operators without the one-body 
operator by extracting the one-pion-exchange (or some other) mechanism, which is nearest 
the one-body operator, from the wave function and connecting it to the one-body operator; a 
two-body operator is formed in this way. The RG running of this operator does not generate 
the kink structure. However, we use the one-body operator here because it is often used in 
practical calculations and the kink part will contribute only marginally to the RG running 
for the reaction near threshold (|qr| ~ 0) which is of our interest here. (At \q\ = 0, the kink 
does not appear.) When fitting the counter terms to the RG low-momentum operator, we 
simply omit the kink structure which obviously cannot be simulated by counter terms. As 
we said, the kink part contributes only marginally to the RG running, and thus we do not 
consider this omission to be influential on a conclusion of this work. 
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IV. RENORMALIZATION GROUP RUNNING, CHIRAL SYMMETRY AND 
COUNTER TERMS 



We study the RG running, guided by the WRG equation, of the n production operator 
presented in the previous section. In %PT, the running coupling constants of the chiral 
counter terms should capture the WRG running of the operator. The leading order counter 
terms (with one spatial derivative) from the heavy-baryon %PT Lagrangian, relevant to this 



work, are given as 22] 

£ct = l^jff2 {bx a N\r • tt<7. V V r • n)NN^N 

+D lh N ] T ■ -kaN ■ JV f (V - V)N + D 1c N\t • tt V - V t • -k)N ■ N ] aN 
+D ld NW ■ irNN^a ■ (V - V)iV + ib le e abc N\r ■ 7T Va + Va T • 7T") <7 b N <J C N 
+iD lf e abc N^T ■ *ka b NN^a c (Va + V a )N 

+bi g N j [r ■ ■ka- V +B- V t ■ tt, r^NN^TiN + D^iV^r • Tra, rJiV • iV f (V + V)r^ 
+?Z)i i e abc Art[ T ■ tt Va - Va r • tt, r^abATArt^^ A^j , (9) 

where D\ a -Du are dimensionless coupling constants. These counter terms are next-to-next- 
to-leading order (N 2 LO) terms in the chiral s-wave pion production operator. 

We explicitly calculate the WRG running of the coefficients of the leading chiral counter 
terms (Di a -Du in Eq. (jHJ)) for an infinitesimal reduction of the cutoff; the LO chiral nuclear 
force is used together. The details of the calculation and the result are given in the App endix. 



This calculation also illustrates that no CSV terms (such as Eq. (25) in Ref. [20(] 6 ) are 
generated in the WRG running. Here, we generalize the result and exclude a possibility 
that the WRG running generates a CSV term, provided that the starting operator does not 
include a CSV term. In Eq. ([3]), the second, third and fourth terms, which we will refer 
to as the generated terms, are to be captured by the counter terms. This part is linear 
in the starting operator. From Eq. ([6]), we see that the starting operator is proportional 
to the pion energy and so are the generated terms. Thus, the counter terms include the 
time derivative of the pion field. Meanwhile, the starting operator and the full Hamiltonian 
have symmetries such as hermiticity, rotational invariance and parity, and the RG running 
maintains these symmetries. To be consistent with the symmetries, each of the generated 
terms has to include a pseudoscalar factor formed by the spatial derivatives and the nucleon 
spin operators. In fact, the heavy-baryon Lagrangian always includes such counter 



6 Discussion in the first paragraph of Sec. IIIII gives a reason why the authors of Ref. 20] obtained the CSV 
term. 
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terms which are the most general terms consisting of the time derivative of the pion field 
and four nucleon fields. Thus no CSV terms are generated in the WRG running. 

Using the result in the appendix, we can (roughly) infer the size of the running in case 
the cutoff is reduced by a certain finite amount. Consider the running of the counter terms 
driven by Owt' [Eq. (J6])] and the one-pion-exchange potential (Vopep) [Eq. (lA.lj) ]. For an 
infinitesimal reduction of the cutoff, the running of the counter terms is given by Eq. (1A.6I) . 
Integrating over A from A down to A (A > A) gives 

In case A > A ^> m^, we have 

Ad-ioop) „ 4^M 3 (A-A) 
la ~ 3(4irf w ) 2 AA ' 1 ; 



The relative strength among Di a and other D\ x (x = b, is the same as Eq. ( 1A.6I) . This 



is the running of the counter terms generated by the one- loop formed by Owt' and Vopep, in 
which the internal relative momentum of the two nucleons runs from A to A. The running 
of the counter terms also captures contributions from the sum of the ladder diagrams in 
which Owt' is dressed by multiple iteration of Vopep- Such a contribution roughly gives the 
running coupling constant as 



A(multi— loops) loop) i ... ... 

U \a ~ U Xa 1 H Ta x- \9 1" 



M(A-A) / M(A — A) 



2 



= Ot"-" (— , (12) 

where A — A < Anf^ is assumed from the first to second line. Although not explicitly shown, 
each term on the R.H.S. of the first line, in fact, is accompanied by a dimensionless constant 
factor which is 0(1) and is different for different terms. Even though Eq. (Tl2|) is a rather 
crude estimation of /)j™ ultl_loo P s ) ? ^ s ^[\\ implies a condition for the WRG running of the 
chiral counter terms to be of natural size: A — A 4-71"/^. (We supposed that the starting 
chiral counter terms are of natural size.) If this condition is not satisfied, the running of 
the counter terms may not necessarily be finite. 7 Furthermore, we find that A — A 47r/ w 
is a condition to justify the naive expectation that the WRG running has the dependence 



7 Our argument here has a similarity to that used for identifying the chiral symmetry breaking scale (A x ) 



with A x ~ 47r/ T 23(. Both of the arguments are based on the size of the RG running. However, we do 



not think that there is a direct connection between the two arguments. In Ref. 23 1, a cutoff variation 



of a reasonable amount leads to a renormalization of couplings, and the size of the renormalization 
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on the cutoff as oc 1/A. For a perturbative calculation, this expectation is true, irrespective 
of the size of A — A, as shown in Eq. In a non-perturbative calculation such as NN 

interaction, however, the scaling of the coupling is more complicated as seen in Eq. (fT2l . 
Now we estimate the size of the generated coupling constant which should be of natural size 
for the counting rule to work. Taking a physically reasonable range (A = 800 MeV and A 
= 500 MeV), which is fairly consistent with the condition obtained above (A — A 4Trf n ), 
we obtain with Eqs. (fTTj) and (fT2l) 

£)(1-1°°P) _ 9) ^(multi-loops) _ ^ ^ (13) 

which are the running of natural size. The above analysis is for the leading order chiral 
counter terms. A similar analysis can also be applied to higher order counter terms to show 
that the running of the couplings is of natural size. Thus this qualitative analysis supports 
the consistency between the RG running and the counting rule in case the cutoff is varied 
within the physically reasonable range. 

In this paragraph, we discuss an implication from the above analysis for the nuclear force. 
Consider the one-pion-exchange potential (Vopep)- The WRG running of Vopep generates 
a series of contact interactions with even numbers of derivatives, i.e., 

VcT(k',k) = ^- 2 (do + ^(e + k") + ..>j , (14) 

where Ci (i = 0, 2, • • •) are dimensionless couplings. These couplings have a RG running 
similar to Eq. ffl2]) . Therefore, if we reduce the cutoff within the reasonable range (e.g., A = 
800 MeV and A = 500 MeV), the resultant renormalization of the couplings is 0(1), which 
implies consistency between the WRG running and Weinberg's counting. However, the 
couplings could be divergent when the cutoff is varied by A — A > 4:nf n . In fact, divergences 
of the RG running of contact AW interactions for some partial waves are indeed observed 
in Ref. j^J (e.g., see Fig. 2 of the paper) in which the cutoff is varied over a very wide range 
from 2 fm _1 to 20 fm _1 ; Eq. ([121 is consistent with the divergent behavior. Obviously, the 
counting rule does not work for the divergent contact interactions. We vary the cutoff within 
the physically reasonable range where the counting rule is not broken by the divergence, and 
study the consistency between the RG running and the counting rule. 

Although the above analysis based on the (semi-) analytic result suggested that the RG 
running of the chiral counter terms is of natural size, it was a rather crude estimation. Thus 

(~ 1 / (4tt /tt) 2 ) is related to the chiral symmetry breaking scale; the RG running is due to one-loop 
diagrams of the irir interaction (perturbative). In our case, we found the size of the cutoff variation 
(~ 4tt/ w ) which may take the couplings away from the natural size; the RG running here is due to the 
sum of the ladder diagrams of the NN interaction (non-perturbative). 
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it is desirable to confirm this result by a numerical calculation as follows. Starting with the 
chiral NLO 7r-production operator [Eq. (jfJJ)] with A = 800 MeV, we reduce the cutoff down to 
A = 500 MeV, with the use of the integral form of the WRG equation [Eq. (j3J)]. 8 Then we 
examine whether the generated terms as a result of the WRG running are captured by the 
chiral counter terms in a way consistent with the counting rule. This time, we use the WRG 
equation [Eq. (J3J)] for which the partial wave decomposition has been done. Because we are 
interested in the AW — > dn reaction near threshold, we only need to consider the 3 Pi — > 3 Si- 
3 Di transition in the AW system. Since we use Eq. (jSJ), the resultant RG effective operator 
does not have spin and isospin structure any more, and only the radial dependence is at 
hand. In fact, some of the chiral counter terms with different spin and isospin structures 
give the same radial dependence in our case ( 3 Pi — > 3 S\- 3 D 1 transition near threshold). 
Thus we know that the generated part includes various counter terms [Eq. Q], but we 
cannot separate them in practice because they give the same radial dependence. What 
we can do is to use a representative of such counter terms to parameterize the generated 
part. Fortunately, this situation still does not ruin our goal of examining the consistency 
between the RG running and the counting rule. We use the following chiral counter terms 
as representatives for parameterizing the generated part: 

£ = Dla p 9A ro UN^t ■ no- V N+H.cAN^N 

+ 2MU\*Mft (iN * T ' ™' ^ N+H - C ^ N ' V N+H.c) 

+ MP {lN]T ' **' * N+H - c -) N]{ y ■ V)(V 2 + vV 



In this work, we consider a pion production near threshold where the on-shell relative momentum is 
about 360 MeV. Therefore, it is impossible for an operator with a cutoff less than 360 MeV to treat the 
pion production. Also, we should consider explicitly degrees of freedom whose details matter to a problem 
under consideration. Therefore, we consider 500 MeV to be an appropriate lowest value for the cutoff. 
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-+2 

-+2 

+ imASWi^ ' 5 JV+ff ' C '»( iV, ( * + V )(i 3t + ^)N + H.e.) 

+ 2i^A^il^ (iiVt ^ * ^ ^ ^+^-^-)(iVt(v - V)(i ^ +^)iV + ^.c.) , (15) 

where the D x (x = la, • • ■ , 5/) are dimensionless coupling constants. This is a power ex- 
pansion with respect to the RG scale A, i.e., in powers of pj A. We presented higher order 
chiral counter terms with three or five derivatives. Although there are also many spin-isospin 
structures possible for the counter terms with three or five derivatives, we presented only the 
representatives which we indeed use in parameterizing the generated part. We do not have 
a criterion to select the set of operators in Eq. (fT5|) . Even if we included some more opera- 
tors, we would not have a better parameterization. The coefficients for those counter terms 
are D 3a -D 3c and D 5a -D 5 f, where the index 3 and 5 indicate how many nucleon derivatives 
they have. Note that contact operators with an even number of spatial derivatives will not 
survive sandwiching between wave functions for the 3 Pi — > 3 Si- 3 Pi transition. 

The counter terms with the LECs L>3 C , D^a, D 5e , and D 5 f are the so-called redundant 
terms. The redundant terms are generated in the WRG running, as explicitly shown in 
the Appendix. Although those terms may be eliminated by a field redefinition, the field 
redefinition is accompanied by a contribution from the measure (the Jacobian factor) which 
is difficult to calculate 24]. Thus, we explicitly consider the redundant terms, which means 



that we explicitly consider the on-shell energy dependence of the low-momentum operator. 



V. NUMERICAL RESULTS 

We start with the chiral NLO s-wave pion production operator for the AW — > dir reaction, 
presented in Sec IIHI The starting operator is defined in the model space with A = 800 MeV 
(sharp cutoff). Using the integral form of the WRG equation in Eq. we calculate the 
RG low-momentum operator for A = 500 MeV. 9 Unless otherwise stated, we use the low- 



momentum AW potential obtained from the CD-Bonn A/W-potential lOj to generate the 
initial proton-neutron 3 Pi scattering and final deuteron wave functions. Some reasons for 
using the phenomenological nuclear force rather than a xPT-based nuclear force have been 
given in the introduction. For the purposes of this paper we can ignore any charge-dependent 



9 We may use either Eq. © after the partial wave decomposition, or Eq. (J3J) to calculate an effective RG 
operator; the result is the same. We use Eq. ([3]) for an easier calculation. 
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FIG. 1: Running of the pion production operator for the 3 Pi — > 3 Si transition in NN — > (for. 
The diagonal momentum space matrix elements are shown. The starting chiral NLO operator 
(A = 800 MeV) is shown by the solid line. After the RG running, we obtain the low-momentum 
operator with A = 700 MeV (dashed line), A = 600 MeV (dotted line), and A = 500 MeV (dash- 
dotted line); r] = 0.1. The abrupt drop in the curves close to the cutoff is a reflection of our 
treatment of the one-body operator (the kink structure). 

effects and assume that our calculation (of np — ► (for ) applies equally well to (Coulomb- 
corrected) pp — > dir + . Thus we will consistently write NN — > (for. We also employ near 
threshold kinematics, i.e., r] ~ 0.1, where rj = q/m n is the emitted pion momentum divided 
by the pion mass. The chiral counter terms that we use to simulate the generated part of 
the RG low- momentum operator have been presented in Eq. (|15|) . 

A. The 3 Pi — ► 3 S\ transition 

The running of the radial part of the 3 P\ — > 3 Si transition operator (diagonal matrix 
elements) is shown in Fig. [TJ The solid line is the starting NLO chiral operator (before RG 
running). After the RG running, we obtain the RG low- momentum operator shown by the 
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dashed line (A = 700 MeV), dotted line (A = 600 MeV) and dash-dotted line (A = 500 MeV). 
Only the result for the diagonal components is given here — a similar trend is found for the 
off-diagonal components. As stated earlier, we expect to observe a kink structure, which is 
similar to the "jump-up" structure in Ref. [jjj, in the low-momentum operators. Here it 
actually appears up as an abrupt drop close to the cutoff. 

Now we parameterize this low-momentum operator using the NLO operator plus the 
counter terms in Eq. ffl5|) . We omit the kink part when fitting the counter terms. Since 
the low-momentum operator has the on-shell energy dependence, we also need to include 
the redundant terms in our fit. The low- momentum operator is calculated for a range in rj 
between 0.02 and 0.1 with steps of 0.02. We fit the LECs to these low-momentum operators 
using the least squares method and with both diagonal and off-diagonal matrix elements 
as input. The //-dependence is then parameterized by the redundant terms. All momenta 
smaller than A are used with equal weight in the fitting. 10 The LECs obtained in this way 
are presented in Table fl] We notice that D\ a changes significantly when adding the D 3x 
terms. This is a consequence of introducing the redundant term D 3c . A similar change is 
observed also in D 3a , D 3 b, and D 3c as the D 5x terms are added. The LECs in Table [J are 
used to obtain the results in the following figures and tables. 

In Fig. [21 we show the simulation of the low-momentum operator using the NLO operator 
plus the counter terms. The dash-dotted line is the starting NLO chiral operator and the 
solid line is the low-momentum operator. The dashed line contains the NLO operator plus 
the counter term with coefficient D\ a . The dotted line contains the additional terms D 3a , 
D 3 b, and D 3c of Eq. (fT5l) . The generated part (the difference between the solid and dash- 
dotted lines) is thus accurately captured by the higher order counter terms. In addition, 
the size of the dimensionless couplings is natural as can be gleaned from Table [J where we 
also include five derivatives counter terms. 11 The natural sizes of the LECs show that the 
chiral expansion of the contact operator indeed converges very well, confirming the visual 
information in Fig. [2j 

Because of the large momentum mismatch between the initial and final on-shell momenta 

10 In principle the smaller momentum should have larger weight in the fitting, since in an effective theory 
the low-momentum behavior of the operator is more accurately described. However, as we will see shortly, 
we can easily achieve a good fit for the full range < A. Therefore we do not think it necessary to put 
heavier weight on the smaller momenta. 

11 The LECs D$ c and D^f are relatively large, although not disturbingly so. However, because of the 
operator structure, the D^c-term appears with a factor 1/3 in the 3 Pi — > 3 Si transition matrix element. 
The £>5/-term has a factor A/M, compared to the IJ^-terms (x = 5a, 56, 5c), which makes its contribution 
smaller by a factor of ~ 2. Therefore, the seemingly large LECs are ascribable to the result of the 
coefficients which accompany the LECs. 
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FIG. 2: Simulation of the low-momentum pion production operator (A = 500 MeV, r\ = 0.1) for the 
3 Pi — » 3 Si transition in NN — » dir. The diagonal matrix elements are shown. The low-momentum 
operator (solid line) is simulated by the original NLO operator (one-body and rescattering terms, 
dash-dotted line) plus a contact term with one nucleon derivative (dashed line). The dotted line 
has additional contact terms with three derivatives. 

in NN — > d,TT, the off-diagonal components are actually more important than the diagonal 
ones. The results for the off-diagonal matrix elements are shown in Fig. [31 The low- 
momentum RG operator is again well parameterized using the higher order counter terms. 
A similar trend is seen also for the other off-diagonal components. 

B. The 3 P 1 -> 3 D 1 transition 

We also examine the running of the operator for the 3 Pi — > 3 P>i transition. The result is 
shown in Fig. HI where the diagonal matrix elements of the operator are given. A counter 
term needs to have at least three derivatives for this transition, so the counter term (Di a ) in 
Eq. (I15p does not contribute. We simulate the low-momentum operator using the counter 
term {Dsb) — the D 3a and D 3c terms cannot contribute. We choose the LECs in the same 
way as described in the previous subsection. This fit can be done independently of the fit of 
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TABLE I: The coupling constants D x for A = 500 MeV. The first column is the transition which 
the counter terms induce. In the second column, "1" means that we use the counter term with one 
nucleon derivative in simulating the low-momentum operator. For "3" ("5"), we additionally use 
the counter terms with three (three and five) derivatives. 

#ofV Dig D 3 q t) 3b D 3c D 5a 

1 0.77 - 

3 Pi 3 5i 3 0.39 -0.64 0.70 1.31 

5 0.42 -0.24 1.98 0.44 0.52 
3 - - -1.37 
5 - - -1.21 

#ofV D 5b D 5c D 5d D 5e D 5f 

3 Pi -» 3 5i 5 5.63 7.63 2.39 -2.08 -8.19 

3 Pi^ 3 L»i 5 0.93 -1.62 - - -1.58 



the 3 Pi — > 3 Si transition operator since the other counter terms not shown in Eq. ffT5| now 
contribute in a different linear combination. Already the one-parameter-fit (dashed line) 
is fairly accurate. When we include the counter terms with five derivatives, we obtain an 
almost perfect simulation as shown by the dotted line in Fig. HI Also these D x are natural 
(see Tabled]). The same level of accuracy is also observed in the simulation of the off-diagonal 
components. 

C. Matrix Elements 

We present in Table ILT1 the matrix elements of the parameterized low-momentum operator 
to show the accuracy of the parametrization using the counter terms. It is clear that the 
counter terms accurately absorb the generated part of the low-momentum operator — already 
with the first counter term the expansion deviates from the low-momentum operator with 
only 1%. With further higher order counter terms included, the matrix elements of the 
low-momentum operator and its parametrization become practically indistinguishable. The 
convergence of the counter term expansion is thus very good. 

In Table (Till the values for the matrix element of each component of the operator are given 
to show the significance of the contribution from the generated terms. At A = 500 MeV, 
the generated part contributes by as much as 11%. This significant contribution for the 
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FIG. 3: Simulation of the low-momentum pion production operator (A = 500 MeV, 77 = 0.1) for 
the 3 Pi — > 3 Si transition in NN — ► dir. The off-diagonal momentum space matrix elements are 
shown for A;fix = 10 MeV. The diagonal components are also shown by the lower solid line. The 
other features are the same as in Fig. [2l 

relatively small cutoff (A = 500 MeV) can be expected, because the cutoff is fairly close to 
the on-shell momentum of the initial NN state (p ~ ^Mm n m 360 MeV). Thus, NN states 
which considerably contribute to the reaction have been integrated out and their strength 
shifted to contact terms. 



D. Wave-function dependence 



In this subsection we investigate the sensitivity of our result to variations in the AW 



potential. Using the AV18 



and the Nijmegen I 26] potentials, instead of the CD-Bonn 
potential, we again study the running of the chiral NLO operator. The result is shown in 
Fig. [5] together with the earlier result for the CD-Bonn potential for comparison. Obviously, 
the NLO operator evolves to different low-momentum operators, depending on the choice of 
nuclear potentials. This result is understandable since the matrix element of the chiral NLO 
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FIG. 4: Simulation of the low-momentum pion production operator (A = 500 MeV, 77 = 0.1) for 
the 3 .Pi — ► 3 Di transition in NN — > dir. The diagonal momentum space matrix elements are 
shown. The low-momentum operator (lower solid line) is simulated by the original NLO operator 
(one-body and rescattering terms, dash-dotted line) plus a contact term with the D^b coefficient 
(dashed line). The dotted line, which almost falls on the solid line, additionally includes counter 
terms with five derivatives. The operator for the 3 Pi — > 3 Si transition is shown for comparison 
(upper solid line). 

operator (A = 800 MeV) rather depends on the choice of the nuclear force, as demonstrated 
in Table IIVI However, the low-momentum operators can be accurately parameterized by 
the chiral NLO operator plus several counter terms with the natural strength, just as we 
have seen in the case of the CD-Bonn potential. The accuracy of the parametrization is at 
a similar level regardless of potentials. 

The discrepancy between the matrix elements calculated by different potentials can be 
regarded as a higher order effect. This is supported by the size of the differences in the 
cross section (Table HV|) . and is also confirmed by the operators — the spread in the curves of 
Fig. [1] stemming from RG running is similar to the spread in Fig. [5] between different wave 
functions. 
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TABLE II: Matrix elements of the pion production operator including the one-body, rescattering 
and several counter terms (77 = 0.1). In the first row, "0" includes no counter terms, i.e., it refers to 
the original chiral NLO operator. "1" includes the counter term with one nucleon derivative, while 
"3" ("5") includes additional counter terms with three (three and five) derivatives. All matrix 
elements are divided by the matrix element of the RG low-momentum operator with the kink 
omitted, to which the counter terms are fitted. 



A (MeV) 





1 


3 


5 


700 


1.02 


1.01 


1.00 


1.00 


600 


1.05 


1.01 


1.00 


1.01 


500 


1.11 


1.01 


1.00 


1.00 



TABLE III: Relative contributions to the matrix element from each component of the pion pro- 
duction operator for the original operator (A = 800 MeV) and a few other cutoffs; 77 = 0.1. The 
symbols, "IB", "WT"' and "Generated term" denote contributions from the one-body, modified 
WT term, and generated terms, respectively. The sum of all contributions is normalized to unity. 



A (MeV) 


(IB) 


(WT') 


(Generated term) 


800 


0.07 


0.93 





700 


0.07 


0.95 


-0.02 


600 


0.05 


1.01 


-0.05 


500 


0.01 


1.09 


-0.11 



In the numerical analysis presented above, the counter terms allowed by the chiral sym- 
metry accurately simulate the generated part of the RG low-momentum operator, with the 
coefficients of the natural strength. We have examined the accuracy of the fitting both 
graphically and quantitatively (Table ILT|) . Thus, our result allows us to conclude that the 
RG running is consistent with the counting rule for the cutoff changed over the physically 
reasonable range. This result is also consistent with our (rough) semi-analytic estimation 
presented in Sec. [TV] 

VI. SUMMARY 

In nuclear %PT with the cutoff regularization, the cutoff has a physical meaning and 
its choice is not arbitrary, however, it still may be varied within a reasonable range. We 
examined whether the running of the counter terms, which is a result of the cutoff variation, 
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FIG. 5: Dependence of the running of the pion production operator ( 3 Pi — > 3 Si) on the choice 
of nuclear force for diagonal momentum space matrix elements. The original chiral NLO operator 
(solid line) evolves to the dashed (AV18), dash-dotted line (Nijmegen I), or the dotted (CD-Bonn) 
line. The cutoff value for the low-momentum operators is A = 500 MeV; rj = 0.1. 



TABLE IV: The reduced s-wave cross section, a = a/r]\ v= o, for various NN potentials. 

CD-Bonn AV18 Nijmegen I 
a (fib) 237 205 224 

is consistent with the counting rule. This consistency is a condition for %PT to be useful. We 
proposed to use Wilsonian RG (WRG) equation for this investigation. With the use of WRG 
equation, it is guaranteed that no chiral-symmetry-violating (CSV) terms are generated, 
provided that we start with an operator consistent with the chiral symmetry. Besides, the 
RG running of the operator, which exactly reflects the high momentum states integrated 
out, is at hand, and thus we can study the convergence of the chiral expansion of the RG 
running. 
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As a demonstration, we applied the WRG equation to the s-wave 7r-production operator. 
We started with the chiral NLO operator 11| with the sharp cutoff, A = 800 MeV. We 
showed that no CSV terms are generated in the WRG running by analytically calculating 
the WRG running of the pion production operator for the infinitesimal reduction of the 
cutoff; the chiral LO nuclear force is used together. In this calculation, we used differential 
form of the WRG equation without the partial wave decomposition. We made the rough 
estimate based on the result of this analytic calculation, and found a range of the cutoff 
variation for which the WRG running is of natural size. We also found that the WRG 
running may be divergent when the cutoff is changed more than 1 GeV; the WRG running 
does not necessarily have a simple oc 1/A behavior when the cutoff is changed by > 1 GeV. 
Therefore, for the variation of the cutoff over the physically reasonable range, the RG running 
is consistent with the counting rule. This argument based on the analytic calculation is also 
applicable to the nuclear force and other transition operators. Thus we made a remark that 
the WRG running of the nuclear force is also consistent with Weinberg's counting rule, as 
long as we are concerned with a variation of the cutoff within the reasonable range. 

In order to confirm the rather rough estimate given above, we numerically calculated the 
running of the chiral counter terms for the variation of the cutoff over the physically reason- 
able range. The CD-Bonn potential is used in this calculation. The couplings of the counter 
terms were fitted to the RG running of the operator given by the WRG equation. We found 
that the parametrization in terms of the counter terms is indeed accurate. Thus, the matrix 
elements for the low- momentum operator can be reproduced within 1% already by includ- 
ing the NLO operator plus just the first (N 2 LO) counter term, with further improvements 
when higher order counter terms are included. All the couplings of the counter terms have 
natural strength, and thus the running can be safely considered to be of higher order than 
the NLO. We also used different phenomenological NN potentials, and found that the result 
is essentially the same as above. As mentioned in the introduction, this result implies that 
the phenomenological NN potentials used do not contain a significant CSV component, 
and thus would serve as a support of the validity of the hybrid approach. The remaining 
difference in the matrix element of the NLO operator among these NN potentials can most 
likely be captured by the N 2 LO counter terms, some of which would be fitted directly to 
data. Thus our result indicates that, if the cutoff is varied within the physically reasonable 
range, the WRG running and the counting rule are consistent, at least for this reaction and 
to the order we are working. 
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APPENDIX: WILSONIAN RENORMALIZATION GROUP RUNNING OF CHI- 
RAL COUNTER TERMS 

Here we show analytically that the running, based on the Wilsonian renormalization 
group (WRG) equation, of a pion production operator is captured by the renormalization of 
the chiral counter terms, provided that we start with a pion production operator consistent 
with chiral symmetry; the operator vanishes at threshold and the chiral limit. We use the 
differential form of the WRG equation [Eq. (j2J)] to study what operators are generated as 
a result of the RG running. We use here the WRG equation without partial wave decom- 
position in order to easily identify the operator structure of the resultant operators. We 
examine the running of the chiral NLO s-wave 7r production operator given in Eq. ([HD- For 
the NN interaction, we use the chiral LO NN interaction, given as: 

n<* , m = . T2 1| '^f + cm + cm 

(A.l) 

with k[ — ki = &2 — k' 2 = k' — k; fei(fc^) is the momentum of the i-th nucleon before (after) 
the interaction. The first term is the one-pion-exchange potential (Vopep) an d the second 
and third terms are contact interactions with different dependences on the nuclear spin. 
The couplings of the contact interactions are functions of A. We start with the set of the 
NN potential and the 7r production operator specified above. The operators are defined in 
the model space with the cutoff A. We reduce A by an infinitesimal amount \8A\, and see 
the running of the 7r production operator, 50. In the following, we replace 0\ in Eq. ([2]) 
with either Owt' or Oib, and V\ with either the Vopep or one of the contact interactions 
(we use the C^-term only), and calculate 50 for each combination of O and V. We will 
see that 50 is captured by the renormalization of the chiral counter terms, except for the 
kink part from the running of Oib- We show this only for the leading counter terms which 
have one spatial derivative. We do not explicitly calculate further higher order counter 
terms from 50. Also, our calculation of 50 does not exhaust every possible combination 
of O and V. Nevertheless, our presentation is enough for our purpose of demonstrating 
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how the RG running is captured by the chiral counter terms. As stated in the text, no 
chiral-symmetry-violating terms are generated by the RG running anyway. We perform 
the following calculation for the threshold kinematics. Above the threshold, the starting 
operators in Eq. have the dependence on the center-of-mass momentum even when we 
set the initial state in the center of mass system, and thus the resultant generated term also 
has the dependence on it. This dependence leads to an additional renormalization of the 
chiral counter terms. 



1. O = Owt' and V = Vqpep 



60£\k>,k;p>,p) = -\5A\M(-A) ( 9 -^ 




( <tx ■ (k[ - A)<ti ■ (A - fci) <x 2 ■ (A + fc 2 ) b 
lj< 2 + (k[ - A) 2 ][m 2 + (A - fci) 2 ](l -p 2 /A 2 ) abc 1 21 1 ' 2j 

o-i ■ (fcj - A) q-2 ■ (-fca ~ A)o-i ■ (A - fcQ ^Ufi«oWA9\ 
[m 2 + (fc^ - A) 2 ] [774 + (A - fci) 2 ](l -p' 2 /A 2 ) 1 Z J 

where the factor (— \8A\) means the infinitesimal reduction of the cutoff. 50^ is captured 
by a series of the counter terms. In order to see this more clearly, we expand the factors in 
the denominator as follows: 

1 1 



ml + (k - A) 2 ml + A 2 + k 2 - 2A ■ k 

1 k 2 - 2 A ■ k 



ml + A 2 (m 2 +A 2 ) 2 



+ •••, (A.3) 



1 



1-P 2 /A 2 = 1+ P 2 /A 2 --- - > ( A - 4 ) 

and so on. The expansion in Eq. (1A.4j) generates the on-shell momentum dependent terms 
which are to be captured by the redundant counter terms. After the expansion of Eqs. (1A.3j) 



and (IA.4I) . we perform the angular integral of A, eliminating the terms with odd numbers 
of A. Here, we just keep terms contributing to the renormalization of the leading counter 
terms shown in Eq. 0. That is, we retain the leading terms of the expansion in terms of 
l/(m 2 + A 2 ) in Eq. (1A.2[) . and take the first term in Eq. (IA.4I) . Note that the second term 
in the second line of Eq. (IA.3I) also gives a contribution. Thus we have, 

_a \5A\M ( g\ \ fg A m n " 



80^(k',k;p',p) 



67r 2 (m 7r 2 + A 2 ) V 4/ 2 y/ V 4/3 J 
{-2a x ■ {k x - k[) + o"! • (feg - k' 2 )}e abc r^r. 



2 



+ {2 CTl ■ (kt + fci) - <t 2 ■ (fci + fe' 1 )}(2i)(r 1 a 



To") 



+ (!«-> 2). (A.5) 
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The terms generated in this way are captured by the renormalization of the chiral counter 
terms in Eq. as follows: 

5D la = 25D lb = -25D lc = -5D ld = 26D lg = -A5D lh = - -Mj^j . (A.6) 

3 (47r/ 7r ) 2 (m^ + A 2 ) 

2. = WT / and V = C s 

} (2)/,/ / \ , XA | hjrr , / \\ ( 9A?n n \ f dVl k ( ax ■ (k[ - A) 



60 A 2 \k',k; P ', P ) = -\SA\MC S (A) ^ / 



4/3 ; 7 (2vr)3 VK 2 + (fei - A) 2 ](l -p 2 /A 2 ) 

o-i • (A - 

K 2 + (A-fc 1 ) 2 ](l-p' 2 /A 2 ) i 



+ r-,2 ■ /r^y I^Hn l + (1 - 2) . (A.7) 



As done for SO A \ we keep the leading contributions [l/(m^ + A 2 )] to obtain 

SOZHk'M = ffgff^ (^) ■ (fa - W^if + (1 « 2) . (A.8) 

This generated term is captured by the renormalization of the chiral counter terms: 

~ C S (A)M*\5A\ 
Wl ° = 247r»« + A») • (A ' 9) 

3. O = Oib and V = Vopep 

50 A 3) {k',k;p',p) = -|5A|M 



#1 \ f ~i9Am 7T (2n) 



Aft) V 4M/ 7r J 
dQ - ' *»> (k' 2 + A) ^ • ^ + • ( A ~ fcl) ^ (A + fc2) Tf ( Tl ■ r 2 ) 



(2tt) 3 V [m 2 + (A-fc!) 2 ](l-p 2 /A 

+ * ,3, <- a - « "'^asy^ H + (i - 2),a - io) 

We perform the angular integral (A) to obtain 

50 A (k,k, P , P )- A2 ^ Ap J ^ j 

X 1^ 2 " A) m 2 + (fci + fc l} 2 ri(TrT2) 

r/ a , v -( fc/ l + fc 2 )Q- 2 -(-fc^ + fc 2 )cr r (-fc 2 + fc 1 ) \ ^nWAin 

+ 5(A - A*) m 2 + (fc , + k2? fa-^ J+(l~2) (A.ll) 

We obtained the operators with the 5-function. This term gives the kink structure found 
in the text, and is not captured by the chiral counter terms. We will see later that the RG 

(3) 

running of 50 A generates the chiral counter terms. 
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4. = Oib and V = C S 
S0^ M = _|, A |M 0i (A)(=^!)/^(^ + A)£^^ 

+i ., ( _ A _ fe!) ^±M) Tr + (1 „ 2) 



|5A|MC 5 (A) (-ig A m. 



{S(k' 2 — A)<Ti • (fc^ — fc 2 



A 2 V *M/ W 

+<5(A - fc 2 ) o-! • (fcx - fc 2 )) 7? + (1 ~ 2) . (A. 12) 



5. O = 5C»i 3) and F = C 5 

In Eqs. (I A. lip and ( 1A.12I) . we obtained the operators with the S- function. The use of the 
WRG equation [Eq. (j2|)] once more gives the counter terms as follows: 

SO f(k ' M - -M Cs( A),-)^ ) f£==) / $ 

(J ' ^ 2 ~~ ^ K + (fe , +A) 2 ](1 _ p2/A2) 

gl -(fej-A) g2 -(-fc 2 -A) < r 1 -(2A) a 
K + (fei-A) 2 ](l-p 2 /A 2 ) 11 2jTl 
o-! ■ (2A)o-i ■ (A - fei) <t 2 ■ (A + fc 2 ) 

[m 2 + (-A + fci) 2 ](l -p ,2 /A 2 ) T ^ Tl ' T2) 

I r a I r/ a , S Q-r(A + fc 2 )<T 2 -(A + fc 2 )o- 1 -(fc 1 -fc 2 ) A 

+ dA d(A - fc 2 J (ri • t 2 JTj + (1 <-> 2) (A.13) 

[m 2 + (A + fc 2 ) 2 ](l -p' 2 /A 2 ) J 

We perform the angular integral and keep the leading terms [oc l/(m 2 + A 2 )] to obtain 

x '\5A\5(k' 2 - A){a x ■ (k[ - k' 2 )a x ■ <t 2 }t?(t x ■ r 2 ) 

+ \5A\8{A - k 2 ){(Ti ■ <r 2 * x ■ (kt - k 2 )}(r x ■ r 2 )rf 

- {2<r 2 • (fei - fei) + 6<r 2 ■ (fc 2 - k' 2 ) - 2i(ffi x <r 2 ) ■ (fei + fe' 1 )}i(r 1 x r 2 ) a 

+ {2<r 2 • + fei) + 6<r 2 • (k 2 + fc' 2 ) - 2i(<n x tr 2 ) ■ - fci)}r 2 a ] + (1 <-> 2) (A.14) 

The terms without 5-function contribute to the renormalization of the counter terms as 
follows: 

8D la = 35D xb = 35D xf = -25D lg = 65D lh = -65 D u = l^^M^M . ( A .l5) 

87r 2 (m 2 +A 2 ) 
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